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ABSTRACT 

We present a comparison of azimuthally averaged radial surface brightness ji(r) 
profiles and analytical bulge-disc decompositions (de Vaucouleurs, r 1 / 4 bulge plus 
exponential disc) for spiral galaxies using Hubble Space Telescope/ Advanced Camera 
for Surveys V-band imaging from the Space Telescope A901/2 Galaxy Evolution Sur- 
vey (STAGES). In the established classification scheme, antitruncated /i(r) profiles 
(Type III) have a broken exponential disc with a shallower region beyond the break 
radius r^k- The excess light at large radii (r > r^k) can either be caused by an outer 
exponential disc (Type Ill-d) or an extended spheroidal component (Type III-s). Using 
our comparisons, we determine the contribution of bulge light at r > rbrk for a large 
sample of 78 (barred/unbarred, Sa-Sd) spiral galaxies with outer disc antitruncations 
(^brk > 24 mag arcsec -2 ). In the majority of cases (~ 85 per cent), evidence indicates 
that excess light at r > rbrk is related to an outer shallow disc (Type Ill-d). Here, the 
contribution of bulge light at r > rbrk is either negligible (~ 70 per cent) or too little 
to explain the antitruncation (~ 15 per cent). However in the latter cases, bulge light 
can affect the measured disc properties (e.g. /ibrk, outer scalelength). In the remaining 
cases (~ 15 per cent), light at r > rbrk is dominated by the bulge (Type III-s). Here, 
for most cases the bulge profile dominates at all radii and only occasionally (3 galaxies, 
~ 5 per cent) extends beyond that of a dominant disc and explains the excess light 
at r > rbrk- We thus conclude that in the vast majority of cases antitruncated outer 
discs cannot be explained by bulge light and thus remain a pure disc phenomenon. 
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INTRODUCTION 

The light profiles of spiral galaxies consist of two prin- 
cipal components: an inner, bulge-dominated component; 
and an outer exponentially declinin g disc with some mi- 
nor deviations related to spiral ar ms (de Vaucouleursll l959; 
iFreemanl |l970). However, since Ivan der Kruitl (|l979f T we 
have known that this 'classical' picture fails for most spi- 
ral galaxies, particularly at the faint surface brightness \i 
of the outer stellar disc. We now know that most disc pro- 
files are best described by a two slope model (broken ex- 
ponential), characterised by an inner and outer exponential 
scalelengt h separated by a r elatively well defined break ra- 
dius rbrk (jPohlen et al.l l20o3 ). Many studies have now re- 
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ported (mainly using surface photometry) the existence of 
broken exponen tial discs, or trunca t ions, i n spiral galaxies in 
both the local (jPohlen et al.1 l200ll2007l: IPohlen fe Truiillo 



2006; iBakos et al.1 120081: lErwin et al.1 120081 : iGutierrez et al 
20 111 : iMaltbv et ai1l201lh and d istant z < 1 Universe (jPerez 



20041 : iTruiillo fe Pohlenl 120051 : lAzzollini et all l2008h . Bro- 
ken exponential discs have also been reporte d using re- 
solved star counts on som e nearby galaxies (|lbata et al.l 
120051 : iFerguson et~alll2007l ) . 

These studies have resulted in a comprehensive clas- 
sification scheme for disc galaxies based on break features 
in t he outer disc componen t of their radial 11 profil es (see 
e.g. lErwin et al.ll2005l . 120081 : IPohlen fe Truiillolbood ). This 
scheme consists of three broad profile types (Type I, II and 
III): Type I (no break) - a single e xponential disc extend- 
ing out to several scalelengths (e.g. Bla nd-Hawthorn et al.l 
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2005); Type II (down-bending break, truncation) - a bro- 
ken exp onential disc with a s hallow inner and st eeper outer 
region (jvan der Kruitl Il979l : IPohlen et ai1l2002h : Type III 
(up-bending break, antitruncation) - a broken exponential 
disc with a shallow er region beyond the break radius rbrk 
(jErwin et alJl2005h . 

In the classical picture (simple bulge and disc), the de 
Vaucouleurs, r 1//4 bulge profile dominates in the centre while 
the exponential disc dominates at larger radii. However, the- 
oretically the r 1//4 profile would domin ate again at some low 
surface brightness. lErwin et al.l (|2005l ) suggest that in some 
Type III profiles (up-bending breaks), the excess light be- 
yond the break radius rbrk could be attributed to light from 
the spheroidal bulge or halo extending beyond the end of the 
disc. Consequently, Type III profiles can be separated into 
two distinct sub-classes depending on whether the outer pro- 
file r > rbrk is dominated by a disc (Typ e IH-d ) or spheroidal 
component (Type III-s). lErwin et al.l ([2005) also propose 
that antitruncations with a smooth gradual transition and 
outer isophotes that are progressively rounder than that of 
the main disc, suggest an inclined disc embedded in a more 
spheroidal outer region such as an extended bulge or halo 
(i.e. Type III-s). U s ing this 'ellipse' method, pr evious works 
(jErwin et al.ll2005l . 120081 : iGutierrez et al.ll201lh have found 
that ~ 40 per cent of their Type III profiles are Type III-s. 

However, the ellipse method is limited for face-on discs 
and cases where the outer/inner disc may have different 
orientations and axis ratios . In these insta nces, bulge-disc 
(B-D) decomposition (e.g. lAllen et"al1l2006h provides a use- 
ful tool to determine the contribution of the two major struc- 
tural components (bulge and disc) to the galaxy's light dis- 
tribution and should provide more conclusive evidence. The 
aim of this work is to use B-D decomposition on a large sam- 
ple of 78 outer disc antitruncations from the Space Telescope 
A901/2 Galaxy Evolution Survey (STAGES; iGrav et al.l 
2009) and assess the fraction of Type III profiles that show 
evidence for the excess light at large radii being caused 
or affected by the spheroidal component. This work builds 
on previous studies by using an improved method for the 
classification of Type III-s/III-d profiles (especially for face- 
on discs) and by using a larger more representative sample 
spanning the range of spiral morphologies. 

Throughout this paper, we adopt a cosmology of Ho = 
70kms _1 Mpc _1 , Q A = 0.7, and Q m = 0.3, and use AB 
magnitudes unless stated otherwise. 



2 DATA AND SAMPLE SELECTION 

This work is entirely based on the STAGES data published 
bv IGrav et~aH (|2Q09h . STAGES is an extensive multiwave- 
length survey targeting the Abell 901/2 multicluster system 
(z ~ 0.167) and covering a wide range of galaxy environ- 
ments. Hubble Space Telescope (HST)/ Advanced Camera for 
Surveys (ACS) F-band (F606W) imaging covering the full 
0.5° x 0.5° (rsj 5 x 5Mpc 2 ) span of the multicluster system 
is complemented by ex tensive obs e rvatio ns including pho- 
tometric redshifts from IWolf et~ai] §003). All imaging and 
data are publicly availabl^E In addition to this, all galax- 

1 http:/ /www. nottingham.ac.uk/astronomy/stages 



ies with apparent R ve g& < 23.5 mag and z prio t < 0.4 (5090 
galaxies) were visually classified by seven members of the 
STAGES team into the Hubble types (E, SO, Sa, Sb, Sc, Sd, 
Irr) and their intermediate classes (Gr ay et al. in prep.). 

Our galaxy sample is drawn from lMaltbv et aD (|201lh . 
This consists of a large, mass- limited (M* > 1O 9 M0), vi- 
sually classified (Sa-Sdm, from Gray et al. in prep) sample 
of 327 face-on to intermediate inclined (i < 60 degrees) spi- 
ral galaxies from both the field and cluster environments. 
However, we remove two galaxies for which B-D decom- 
position fails (iVtot = 325). The 182 cluster spirals are at 
a redshift of z c \ — 0.167 and the 143 field spirals span 
a red shift range of 0.05 < z p hot < 0.30. iMaltbv et al.l 
(|201lh analysed the ji{r) profiles for these galaxies in or- 
der to identify broken exponentials in the outer stellar disc 
\i > 24 mag arcsec -2 (their criteria for selecting intrinsically 
similar outer breaks) . Three independent assessors agreed on 
a sub-sample of 78 antitruncated (Type III) outer fi profiles 
(Mbrk > 24 mag arcsec -2 ). We shall use both this Type III 
sub-sample and the total sample in this work. 



3 METHODOLOGY 

For each galaxy in our sample, we perform a two-dimensional 
B-D decomposition based on a two component galaxy 
model comprising of a de Vaucouleurs (r 1//4 ) bulge and a 
single exponential disc. Decompositions were carried out 
on the STAGES F-band imaging using the GALFIT code 
(|Peng et al.1 120021 ) and the method of iHovos etafl (|201lh 
adapted to perform two component fits. Several measurable 
properties are produced for each galaxy including: position 
[x,y], effective radii, total magnitudes, axis ratios, position 
angles for the bulge and disc components, and a sky-level 
estimation. 

B-D decomposition can be sensitive to the initial con- 
ditions used to search the B-D parameter space (e.g. initial 
guess for bulge-to-disc ratio B/D). Therefore, we perform 
two runs of the B-D decomposition with different initial 
conditions from the two extremes. One run starting from 
a bulge-dominated system (B/D = 9), and the other run 
starting from a disc-dominated system (B/D = 1/9). Com- 
parison of these runs (hereafter, Run 1 and Run 2) allows for 
an assessment of the uniqueness/stability of B-D decompo- 
sition on a galaxy-galaxy basis. In the vast majority of cases 
(~ 85 per cent) the results were effectively the same, ~ 70 
per cent being exactly the same and ~ 15 per cent having 
only minor differences that do not affect our analysis (B/D 
the same within ~ 10 per cent). In a few cases (~ 10 per 
cent) the decomposition was catastrophically unstable with 
Run 1/2 yielding both bulge- and disc-dominated systems. 
The remaining cases (~ 5 per cent) showed moderate insta- 
bilities great enough to affect the assessment of bulge light in 
the outer regions of the galaxy. These stability fractions are 
the same in our Type III sub-sample. The unstable solutions 
are mainly driven by differences in the sky level determined 
during the decomposition. However, the overall conclusions 
of this work are not affected by these unstable solutions. 

For each galaxy, we also use the IRAF task ellipse (STS- 
DAS package - version 2.12.2) in order to obtain azimuthally 
averaged rad ial ii(r) profiles f rom the STAGES V-band 
imaging, see lMaltbv et al.l (|201 lh for full details of the ellipse 
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fitting method used. For all our ellipse fits the galaxy cen- 
tre is fixed (all isophotes have a common centre) using the 
galaxy centre deter mined during B-D decomposition. The 
bad pixel masks of Gra y et al . (2009) are also used to re- 
move everything not associated with the galaxy itself from 
the isophotal fit (e.g. light from companion galaxies). 

Using a similar procedure t o previous works 

(jPohlen fe Truiillol I200& lErwin et all l2008h . we fit two 
different sets of ellipses to each galaxy image. The first is a 
free parameter fit (fixed centre, free ellipticity e and position 
angle PA) and tends to follow morphological features such 
as bars and spiral arms. Consequently, these free- fits are 
not suitable for the characterisation of the underlying outer 
disc, however they may be used to determine the e and 
PA of the outer disc component, see iMaltbv et alJ (| 20 111 ) 
for further details. A fixed-parameter fit (fixed centre, e 
and PA using the e and PA determined for the outer disc) 
is then used to produce our final measured ji(r) profiles. 
The necessary sky subtraction is then performed using the 
sky-level estimates generated during B-D decomposition. 
Please note t hat these sky-values sometimes differ slightly 
from those of IMaltbv et~aU (|201lh . 

Analogous fixed-parameter fits (using the e and PA de- 
termined for the outer disc) are also carried out on the disc- 
residual images (ACS image - bulge-only model) resulting 
in a measured \i profile for the disc component /iDisc(^). 
We also obtain azimuthally averaged radial \i profiles for 
the decomposed B-D model using the same fixed-parameter 
ellipses (isophotes) as in the other profiles. These yield sepa- 
rate radial fi profiles for both the disc and bulge model along 
the semi-major axis of the outer disc. 

All resultant \i profiles are corrected for Galactic fore- 
ground extinction, individual galaxy inclination z, and sur- 
face brightness dimming (the \i profiles of field galaxies, 
0.05 < Zphot < 0.30, are corrected to the redshift of the 
cluster z — 0.167). Full details of the fitting procedure, sub- 
sequent photometric calibration, and an estimation of the 
error in the s k y sub traction (±0.18 counts) can be found in 
IMaltbv et ail |201lh . 



4 RESULTS 

B-D decompositions using a de Vaucouleurs (r 1 / 4 ) bulge plus 
an expone ntial disc can be classified into 4 distinct profile 
types fe.g. lAllen 6^11120061 ). see Fig. [2 

(1) Type A: 'Classical' system. The bulge profile domi- 
nates at the centre, while the disc profile dominates at large 
radii. The bulge/disc profiles cross only once. 

(2) Type B: Disc-dominated system. The disc profile dom- 
inates at all radii, with a weak contribution from the bulge 
in the centre. The bulge/disc profiles never cross. 

(3) Type C: The bulge profile dominates at small/large 
radii, but the disc profile dominates at intermediate radii. 

(4) Type D: Bulge-dominated system. The bulge profile 
dominates at all radii with a weak underlying disc compo- 
nent. The bulge/disc profiles never crosfl. 



2 Note the use of capital letters in profile types to avoid confusion 
with other classification schemes. 
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Figure 1. Bulge-disc (B-D) profile types. Top left: Type A, 'clas- 
sical' profile. Top right: Type B, disc-dominated. Middle left: 
Type C, bulge-dominated at small/large radii but disc-dominated 
at intermediate radii. Middle right: Type D, bulge-dominated. 
Bottom left: Type E, 'constrained' outer bulge. Bottom right: 
profile type distributions for the total sample and Type III sub- 
sample. 



In addition to this, we also observe decompositions 
where the disc profile dominates in the centre, while the 
bulge profile dominates at large radii (hereafter Type E). 
Here, we believe an outer antitruncated disk has incorrectly 
affected the bulge profile fit. In these cases, B-D decomposi- 
tion is not a true representation of the galaxy at large radii 
and these galaxies probably have Type B-like compositions. 
Similar constraints may also occur in some Type D profiles. 

Fig. [1] shows the distribution of B-D profile types for 
both the total sample and the Type III sub-sample. As ex- 
pected, the fraction of Type C/E profiles is greater in the 
Type III sub-sample. For the total sample, several other cor- 
relations between B-D profile type, Hubble-type morphology 
(Sa-Sd) and measured B/D ratio are observed. Type A/B 
profiles are equally probable in all Hubble types but Type 
C/D profiles are more common in earlier Hubble types (Sa- 
Sb). Also as expected, mean/median B/D decreases with 
progressively later Hubble types and increases for the se- 
quence of B-D profile types B-A-C-D (increasing bulge 
dominance) . 

For each galaxy in our Type III sub-sample, we com- 
pared the measured, fixed ellipse /n(r) profile with the model 
H profiles from the B-D model in order to assess the con- 
tribution of bulge light in the outer regions of the galaxy 
(r > rbrk)- Fig. [2] shows some examples for each B-D profile 
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Figure 3. A rare/unusual example of a bulge profile causing an 
antitruncation in a \i profile, (Type C example from Fig. [2] break 
radius r^k)- The bulge (red dotted line), disc (blue line), and 
bulge + disc (black dashed line) profiles from B-D decomposition 
are overplotted on the measured /i profile (red circles). The disc- 
residual /iDisc profile (measured \i profile minus bulge-only model, 
red crosses) shows no antitruncation (for errors, see Fig. [2}. 

type. These comparisons resulted in three possible scenarios. 
Bulge light in the outer profile (r > rbrk) either had: 

(i) little or no contribution (~ 70 per cent): For all 
Type A/B profiles the bulge contributes virtually no light 
at r > rbrk and in some Type C/E profiles the contribution 
is negligible. This can be determined by inspection of the 
measured disc-residual profile /iDisc(r) and assessing if the 
properties of the outer profile/break (rbrk, Mbrk, scalelength) 
have been affected with respect to the sky-subtraction error. 
No Type D profiles are in this category. 

(ii) minor contribution (~ 15 per cent): Approximately 
half of these cases are Type C profiles where the bulge profile 
emerges from the end of the disc and contributes some light 
to the outer regions of the galaxy. The remaining half are 
Type E profiles where the bulge appears to be constrained by 
an outer exponential disc. The amount contributed is enough 
to affect the outer profile causing /ibrk and the outer scale- 
length to be different in the disc-residual profile /iDisc(r). 
However, the antitruncation remains present. 

(hi) major contribution (~ 15 per cent): Here, the bulge 
contributes the majority of the light at r > rbrk- Approx- 
imately half of these cases are Type D profiles where the 
bulge dominates at all radii. For these cases, the de Vau- 
couleurs profile is either interpreted as an antitruncation or 
being constrained by an outer shallow disc. The remaining 
half are all Type C profiles (with one exception which is 
Type E). In these latter cases the antitruncation is not ob- 
served in the disc-residual profile /iDisc(r), see Fig.[3]for such 
an example. 



is important to note that for some of these cases (~ 15 per 
cent), bulge light can affect the properties of the disc profile 
(e.g. /ibrk, outer scalelength). For the remaining Type III 
profiles (~ 15 per cent), the excess light at r > rbrk can be 
attributed to the bulge profile (Type III-s). However, few 
of these latter cases (only 3 galaxies with stable decomposi- 
tions, ~ 5 per cent) exhibit profiles where the bulge profile 
extends beyond a dominant disc (Type C) and causes an 
antitruncation in the ^ profile. 

Previous works, (|Erwin et al.l 120081 ; I Gutierrez et al.l 

1 20 111 ) have used the ellipse method to classify their Type III- 
d/III-s profiles. The results of these works are in good 
agreement with their Type III profiles being ~ 60 per cent 
Type Ill-d and ~ 40 per cent Type III-s. These results ap- 
pear not to be consistent with our results. However, consider 
the case where the contribution of bulge light at r > rbrk is 
too little to explain the antitruncation but great enough to 
affect the properties of the profile. The increased contribu- 
tion of bulge light to the /i(r) profile at rbrk could lead to 
a smoothing of the inflection and a roundening of the outer 
isophotes. Therefore, using the ellipse method would likely 
have resulted in these Type Ill-d profiles being classified as 
Type III-s. Applying this reasoning to our results, we obtain 
fractions of ~ 70/30 per cent for Type III-d/III-s profiles 
respectively. Therefore, the ellipse method potentially could 
lead to genuine disc breaks being classified as Type III-s. 
Our method offers an improved way to determine whether 
an antitruncation is disc or spheroid related. 

However, our method does have some obvious draw- 
backs. In a two-component B-D decomposition, an outer 
antitruncated disc could cause the bulge profile to be con- 
strained and lead to an over- estimation of bulge light in 
the outer regions of the galaxy. This naturally enhances our 
fraction of Type III-s profiles which therefore represents an 
upper limit to the fraction of genuine Type III-s profiles in 
our Type III sample. Therefore, we can conclude that in the 
vast majority of cases Type III profiles are indeed a true disc 
phenomenon. 

Several studies have proposed formation scenarios for 
Type H I profiles, mainly via satellite accretion or minor 
mergers f Penarrubia et al.1l2006l : lYounger et alj2QQ7h . A dis- 
cussion of these scenarios is beyond the scope of this work. 
However, we do make one important comment. Pre vious 
works fe.g. lPohlen fe Truiillol[2006l : lMaltbv et al.ll201lh have 
reported that Type III profiles are more frequent in earlier 
Hubble types. This result is consistent with a minor merger 
scenario for their origin. However, the excess Type III pro- 
files in early-types could easily have been related to a natural 
increase in the number of Type III-s profiles in earlier Hub- 
ble types. Fortunately, we observe exactly the same relations 
using just our genuine Type Ill-d profiles and thus our re- 
sults remain consistent with the minor merger scenario for 
the formation of Type III profiles. 



5 CONCLUSIONS 

Our results suggest that for the majority of Type III profiles 
(~ 85 per cent), the excess light beyond the break radius rbrk 
is related to an outer shallow disc (Type Ill-d). However, it 
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Figure 2. Example B-D decompositions and antitruncated profiles. Rows, top to bottom: decompositions producing Type A, Type B, 
Type C, Type D and Type E profiles. The bulge (red dotted line), disc (blue line), and bulge + disc (black dashed line) profiles from B-D 
decomposition are overplotted on the measured /i profiles (red circles). The i/ii/iii notation after the profile type indicates whether the 
bulge profile has a negligible, minor or major contribution outside the break radius r > rb r k 5 as m Section 0] Errors in the measured \i 
pro files are for an over- and undersubtraction of the sky by ±1<7. The /x C rit /Vlim levels represent +lcr/+3cr above the sky respectively, 
see lMaltbv et al.l (|201lh for full details. The model images show the effective radius isophote for both the bulge (red line) and disc (blue 
dashed line) models. 
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